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Abstract. The NA60 experiment at the CERN SPS has measured muon pairs with unprecedented precision 
in 158A GeV In- In collisions. A strong excess of pairs above the known sources is observed in the whole mass 
region 0.2<M<2.6 GeV. The mass spectrum for M<1 GeV is consistent with a dominant contribution 
from tt + -k~ — > p — > fi + fi~ annihilation. The associated p spectral function shows a strong broadening, 
but essentially no shift in mass. For M>1 GeV, the excess is found to be prompt, not due to enhanced 
charm production, with pronounced differences to Drell-Yan pairs. The slope parameter T e « associated 
with the transverse momentum spectra rises with mass up to the p, followed by a sudden decline above. 
The rise for M<1 GeV is consistent with radial flow of a hadronic emission source. The seeming absence 
of significant flow for M>1 GeV and its relation to parton-hadron duality is discussed in detail, suggesting 
a dominantly partonic emission source in this region. A comparison of the data to the present status of 
theoretical modeling is also contained. The accumulated empirical evidence, including also a Planck-like 
shape of the mass spectra at low pr and the lack of polarization, is consistent with a global interpretation 
of the excess dimuons as thermal radiation. We conclude with first results on u in-medium effects. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 12.38. Mh Quark-gluon plasma - 13.85. Qk Lepton Pairs 



1 Introduction 

Dileptons are particularly attractive to study the hot and 
dense QCD matter formed in high-energy nuclear colli- 
sions. In contrast to hadrons, they directly probe the en- 
tire space-time evolution of the expanding system, escap- 
ing freely without final-state interactions. At low masses 
M<1 GeV (LMR), thermal dilepton production is medi- 
ated by the broad vector meson p (770) in the hadronic 
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phase. Due to its strong coupling to the rnr channel and 
the short life time of only 1.3 fm/c, "in-medium" modi- 
fications of its mass and width close to the QCD phase 
boundary have since long been considered as the prime 
signature for chiral symmetry restoration [UEIG] • At inter- 
mediate masses M>1 GeV (IMR), it has been controver- 
sial up to today whether thermal dileptons are dominantly 
produced in the earlier partonic or in the hadronic phase, 
based here on hadronic processes other than tttt annihi- 
lation. Originally, thermal emission from the early phase 
was considered as a prime probe of deconfinement [HIS]. 
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Fig. 1. Background-subtracted mass spectrum before (dots) Fig. 2. Excess dimuons compared to theoretical predic- 



and after subtraction of the known decay sources (triangles) . 



tions [16], renormalized to the data in the mass interval M<0.9 
GeV. No acceptance correction applied. 



Experimentally, it took more than a decade to mas- 
ter the challenges of very rare signals and enormous com- 
binatorial backgrounds. The first clear signs of an ex- 
cess of dileptons above the known decay sources at SPS 
energies were obtained by CERES [HE] for M<1 GeV, 
NA38/NA50 [| for M>1 GeV and by HELIOS-3 [9] for 
both mass regions (see [10] for a short recent review in- 
cluding the preceding pp era and the theoretical mile- 
stones). The sole existence of an excess gave a strong 
boost to theory, with hundreds of publications. In the 
LMR region, 7T7T annihilation with regeneration and strong 
in-medium modifications of the intermediate p during the 
fireball expansion emerged as the dominant source. How- 
ever, the data quality in terms of statistics and mass res- 
olution remained largely insufficient for a precise assess- 
ment for the in-medium spectral properties of the p. In 
the IMR region, thermal sources or enhanced charm pro- 
duction could account for the excess equally well, but that 
ambiguity could not be resolved, nor could the nature of 
the thermal sources be clarified. 

A big step forward in technology, leading to completely 
new standards of the data quality in this field, has re- 
cently been achieved by NA60, a third-generation exper- 
iment built specifically to follow up the open issues ad- 
dressed above [11]. Initial results on mass and transverse 
momentum spectra of the excess dimuons have already 
been published [Tl"lll"2"I[r3"] , supplemented by recent results 
on acceptance-corrected mass spectra and polarization [14] . 
This paper takes a broader view on the results, discussing 
in some detail the observed M-pT correlations, the con- 
nection to hadron-parton duality, and the present status 
of theoretical modeling. Certain aspects of the pr spectra, 
of centrality dependencies and evidence for lo in-medium 
effects are reported here for the first time. 



2 Mass spectra and the p spectral function 

Fig. [1] shows the centrality-integrated net dimuon mass 
spectrum for 158A GeV In-ln collisions in the LMR re- 
gion. The narrow vector mesons u> and <f> are completely 
resolved; the mass resolution at the ui is 20 MeV. The 
peripheral data can be completely described by the elec- 
tromagnetic decays of neutral mesons [T2l[T3] . This is not 
true for the more central data as plotted in Fig. [TJ due 
to the existence of a strong excess of pairs. The high data 
quality of NA60 allows to isolate this excess with a priori 
unknown characteristics without any fits: the cocktail of 
the decay sources is subtracted from the total data using 
local criteria, which are solely based on the mass distribu- 
tion itself. The p is not subtracted. The excess resulting 
from this difference formation is illustrated in the same fig- 
ure (see [T2l[T3l[T5] for details and error discussion). The 
subtracted data for the rj, u> and <f> themselves are subject 
to the same further steps as the excess data and are used 
later for comparison. 

The common features of the excess mass spectra can 
be recognized in Fig.O A peaked structure is always seen, 
residing on a broad continuum with a yield strongly in- 
creasing with centrality (see Fig. [5] below), but remain- 
ing essentially centered around the nominal p pole [15] . 
Without any acceptance correction and px selection, the 
data can directly be interpreted as the space-time averaged 
spectral function of the p, due to a fortuitous cancellation 
of the mass and px dependence of the acceptance filtering 
by the phase space factors associated with thermal dilep- 
ton emission [15] , The two main theoretical scenarios for 
the in-medium spectral properties of the p, broadening [5] 
and dropping mass [3], are shown for comparison, both 
evaluated for the same fireball evolution [TB] . Since agree- 
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Fig. 3. Excess dimuons for semicentral collisions compared to the theoretical model results by Hees/Rapp [18] (left) and 
Renk/Ruppert et al. [19] (right). No acceptance correction applied. 



ment between modeling and data would imply agreement 
both in shape and yield, the model results are normal- 
ized to the data in the mass interval M< 0.9 GeV, just 
to be independent of the uncertainties of the fireball evo- 
lution. The unmodified p, also shown in Fig. [2] (vacuum 
p), is clearly ruled out. The broadening scenario indeed 
gets close, while the dropping mass scenario in the version 
which described the CERES data reasonably well [2,3,6 
completely fails for the much more precise NA60 data. 
A strong reduction of in-medium VMD as proposed by 
the vector manifestation of chiral symmetry [17) would 
make hadron spectral functions in hot and dense mat- 
ter altogether unobservable, but central aspects of this 
scenario are totally unclear, and quantitative predictions 
which could be confronted with data have not become 
available up to today. 

A comparison of the same excess mass spectrum to 
two more recent theoretical developments, covering now 
both the LMR and the initial part of the IMR region, is 
contained in Fig. [3] In contrast to Fig. [21 the theoretical re- 
sults are not renormalized here, but shown on an absolute 
scale. In the p-like region with 7r + 7r~ — > p — > p + pT as the 
dominant source, Hees/Rapp [TS] use the original many- 
body scenario with a p spectral function strongly broad- 
ened by baryonic interactions [2J, while Renk/Ruppert 's 
results [19] are based on the spectral function of Eletsky 
et al. [2U] where the broadening effects from baryons are 
somewhat weaker; that difference is directly visible in the 
low-mass tails of the theoretical mass spectra. The overall 
agreement between the data and the two theoretical sce- 
narios is quite satisfactory in this region, also in absolute 
terms. 



3 Mass spectra and parton-hadron duality 

Moving up into the IMR region M>1 GeV, 2ir processes 
become negligible, and other hadronic processes like 4tt 



(including vector- axialvector mixing in case of [18) ) and 
partonic processes like quark-antiquark annihilation qq — > 
p + p~ take over. The two theoretical scenarios in Fig. [3J 
also describe this part. However, there is a very interesting 
and instructive difference between them. While the total 
yield of the data for M>1 GeV is described about equally 
well, the fraction of partonic processes relative to the total 
is small in [18] where a first-order phase transition is used, 
and dominant in 19J which uses a cross-over phase tran- 
sition. This feature is often referred to as "parton-hadron 
duality" and formed the basis of the successful descrip- 
tion of the NA50 dimuon enhancement in the IMR region 
in terms of thermal radiation |21) . Here, the individual 
sources were not even specified. 

Caution should, however, be expressed as to the use 
of the term "duality" in this context. Parton-hadron du- 
ality is a statement on dilcpton emission rates, dating 
back to the time-reversed process of hadron production 
in e + e~ collisions. It implies that the emission rates us- 
ing either partonic (pQCD) or hadronic degrees of free- 
dom merge together, i.e. become "dual", if the system 
approaches deconfinement and chiral restoration. The va- 
lidity of duality down to masses of 1 GeV, mainly due 
to vector-axialvector mixing, was first shown by Li and 
Gale [22] (see also [2J). However, experiments measure 
yields, i.e. rates integrated over space-time. Duality in 
the yields is not obvious and becomes questionable, if the 
space-time trajectories are different for genuine partonic 
and hadronic processes. Such a difference automatically 
appears through the elementary assumption that partonic 
processes only act "early", i.e. from Tj until T c , while 
hadronic processes (like nir) only act "late" , i.e. from T c to 
thermal freeze-out Tj . If theoretical scenarios are different 
in their trajectories (both as to partonic and to hadronic 
processes), the integrated total yields will, in general, be 
different. Since the scenarios of [TH] and [TH] in Fig. [3] are 
indeed very different (see above), the seemingly equiva- 
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lent description of the data cannot be traced to duality, 
but must be due to internal parameter choices. 

Explicit insight beyond duality, be it real or fortu- 
itous, can be obtained experimentally in the following way 
(see [5] and in particular [13]). In contrast to real photons, 
virtual photons decaying into lepton pairs are character- 
ized by two variables, mass M and transverse momen- 
tum px- Historically, the interest has largely focused on 
mass because of its rich and often structured information 
content, including now the p spectral function discussed 
above. Transverse momentum, on the other hand, contains 
not only contributions from the spectral function(s), but 
encodes the key properties of the expanding fireball, tem- 
perature and, in particular, transverse (radial) flow. The 
latter causes a blue shift of pr, analogous to the case of 
hadrons. In contrast to hadrons, however, which always 
receive the full asymptotic flow reached at the moment of 
decoupling from the flowing medium, lepton pairs are con- 
tinuously emitted during the evolution, reflecting a space- 
time folding over the temperature-flow history in their fi- 
nal pt spectra. Since flow builds up monotonically during 
the evolution, being small in the early partonic phase (at 
SPS energies, due to the "soft point" in the equation-of- 
state), and increasingly larger in the late hadronic phase, 
the final px spectra keep memory on the time ordering of 
the different dilcpton sources, mirrored in a characteristic 
mass dependence of the pt spectra. We shall come back to 
this point below. 

4 Acceptance-corrected mass and p T spectra 

Quantitative insight into the physical meaning of the ex- 
cess dileptons requires a full correction of the data for 
geometrical acceptance and pair efficiencies of the NA60 
apparatus, including the effects of the trigger system. Re- 
sults from Monte Carlo simulations of the acceptance are 
contained in ,13,24;, showing significant variations and in 
particular a strong decrease at low mass and low px- In 
principle, the correction requires a 4-dimensional grid in 
the space of M-pT-y-cosOcs (where Ocs is the polar angle 
of the muons in the Collins Soper frame). To avoid large 
statistical errors in low-acceptance bins, it is performed 
instead in 2-dimensional M-px space, using the measured 
y and cosQ distributions as an input. The latter are, in 
turn, obtained with acceptance corrections determined in 
an iterative way from MC simulations matched to the data 
in M and px ■ The y-distribution is found to have the same 
rapidity width as dN c h/drj, a y ~ 1.5 [24]. The cosQcs dis- 
tributions for two mass windows of the excess and the u> 
are contained in [14]. Within errors, they are found to be 
uniform, implying the polarization of the excess dimuons 
to be zero, in contrast to Droll- Yan and consistent with 
the expectations for thermal radiation from a randomized 
system. 

The outcome for the two major variables M and px 
is first discussed separately for the LMR and the IMR 
regions. Fig. [4] shows a set of mass spectra for some se- 
lected slices in pr to illustrate the evolution from low to 
high px (a pr-integrated mass spectrum over the whole 



mass region is contained in [11]). The spectra are normal- 
ized to dN c h/dr] in absolute terms, using the same proce- 
dure as described in detail for the <p in [25] and relating 
N part ~ dN c h/drj at 77=2.9 as measured to within 10% by 
the Si pixel telescope. Recent theoretical results on ther- 
mal radiation from three major groups working in the field 
are included for comparison 18,19,23,26,27], calculated 
absolutely (not normalized to the data). Results from a 
fourth one [28] only cover the 2ir region and are not yet 
available in p^-differential form. The general agreement 
between data and model results both as to spectral shapes 
and to absolute yields is most remarkable, supporting the 
term "thermal" used throughout this paper. 

At very low pr, a strong rise towards low masses is 
seen in the data, reflecting the Boltzmann factor, i.e. the 
Plank-like radiation associated with a very broad, nearly 
flat spectral function. Only the Hees/Rapp scenario [15J 
is able to describe this part quantitatively, due to their 
particularly large contribution from baryonic interactions 
to the low-mass tail of the p spectral function (as in [28]). 
This was already mentioned in connection with Fig. [3j but 
is much more clearly visible at low px than without any 
Pt selection. At higher pt, the influence of radial flow in- 
creasingly changes the spectral shapes, and at very high 
Pt all spectra appear p-like. However, sizable differences 
between the different theoretical scenarios also exist in this 
region. For example, Hees/Rapp [26] use a hard-scattering 
p which contributes to fill up the p region beyond the 
Cooper-Frye freeze-out p. They also use an extrapola- 
tion of the Drell-Yan process down to the photon point 
(M— >0). While this contribution is small for the whole 
LMR region without pt selection, the low-mass/high-pT 
part is in their case completely dominated by DY. The 
size of radial flow, the major issue here, also varies be- 
tween the three groups. It is recognizably too low in the 
hydrodynamics scenario [S7J , and maximal for the fireball 
scenario of [T9II23] tailored to the NA60 hadron data. 

Fig. [5] (left) shows the centrality-integrated mj spec- 
tra, where tht— (p T + M 2 ) 1/2 , for four mass windows; the 
<fi is included for comparison. The ordinate is absolutely 
normalized to dN c h/drj as in Fig. 2] Apart from a pe- 
culiar rise at low tut (<0.2 GeV) for the excess spectra 
(not the 4>) which only disappears for very peripheral col- 
lisions 10,13 , all spectra are pure exponentials. The rise 
is outside of any systematic errors as discussed in [13]. 
The relative yield associated with it is about 10-20%, 
roughly independent of mass, which excludes a connection 
to the low pt rise seen in pion pt spectra. The absolute 
yield steeply decreases with mass, reminiscent of Dalitz de- 
cays. However, a consistent physical interpretation is still 
open. The lines in the exponential region are fits to the 
data with the function 1/mTdN/dmT oc exp(—mT /T e ff) , 
where the effective temperature parameter T e g is the in- 
verse slope of the distributions. For the excess data, the 
fits are restricted to the range 0.4<pt<1-8 GeV (roughly 
0.1<mr- M<1.2 GeV) to exclude the increased rise at low 
my. Obviously, the slopes depend on mass. Fig. [5] (right) 
shows a more detailed view into the p-like mass window, 
exploiting the same side-window method as used in con- 
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nection with Fig. [5] below, to determine the px spectra 
separately for the p peak and the underlying continuum. 
All spectra are purely exponential up to the cut-off at 
Pt=3 GeV, without any signs of an upward bend charac- 
teristic for the onset of hard processes. Their slopes are, 
however, quite different (see below). 

The central NA60 results in the IMR region [11] are 
shown in Fig. [5] The use of the Si-vertex tracker allows to 
measure the offset between the muon tracks and the main 
interaction vertex and thereby to disentangle prompt and 
offset dimuons from D decays. The offset distribution is 
found to be perfectly consistent with no charm enhance- 
ment, expressed by a fraction of 1.16±0.16 of the level ex- 
pected from upscaling the NA50 results on the IMR in p- A 
collisions 1 lj . The observed excess is really prompt, with 
an enhancement over Drell-Yan by a factor of 2.4±0.08. 
The excess can now be isolated in the same way as was 
done in the LMR region, subtracting the measured known 
sources, here DY and open charm, from the total data. 
Fig. [6] shows the decomposition of the total into DY, open 
charm and the prompt excess. The mass spectrum of the 
excess is quite similar to the shape of open charm and 
much steeper than DY; this explains of course why NA50 
could describe the excess as enhanced open charm [8]. 

The transverse momentum spectra are also much steep- 
er than DY. Moreover, the spectra depend on mass and 
do not show the factorization between mass and px char- 



acteristic for DY, where a common Gaussian distribution 
with a fixed sigma kx=0-8 GeV describes all px spec- 
tra independent of mass. The transverse mass spectra are 
shown in Fig. [7] for three consecutive mass windows. All 
spectra are essentially exponential. However, the steep- 
ening observed at very low my in the lowest mass win- 
dow, seen already before for all masses in Fig. [5] including 
this window, seems to be switched-off in the upper two 
mass windows. As in Fig. [5] the lines are exponential fits 
to the data, restricted again to px>0-5 GeV to exclude 
the rise at low mx- The extracted inverse slope parame- 
ters are 199±21(stat)±3(syst), 193±16±2 and 171±21±3 
MeV, respectively, i.e. about the same within the (rather 
large) errors. 

The inverse slope parameters T e g extracted from the 
exponential fits to the mx spectra are plotted in the left 
panel of Fig. [8] vs. dimuon mass, unifying the data from 
the LMR and IMR regions. In the LMR part, a finer bin- 
ning is used than in Fig. [5] and the coarser-binned data 
(contained in [TJ]) are left out for clarity. For M<1 GeV, 
a correction for Drell-Yan pairs is not done, due to their 
small contribution [13] and the intrinsic uncertainties at 
low masses In the extended LMR analysis up to 1.4 
GeV, the 2 (square) points are corrected, as are all points 
of the IMR analysis (see above). In the region of overlap, 
the data are not statistically independent. The hadron 
data for 77, lu and (f> obtained as a by-product of the cock- 
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tail subtraction procedure are also included in Fig. [5J as 
is the single value for the p-peak from the right panel of 
Fig. \5\ Interpreting the latter as the freeze-out p without 
in-medium effects, consistent with all present theoretical 
modeling [18„ 19 , 23 , 261127] . all four hadron values together 
with preliminary 7r~ data from NA60 can be subjected to 



a simple blast wave analysis [14] . This results in a reason- 
able set of freeze-out parameters of the fireball evolution 
and suggests the following consistent interpretation for the 
hadron and dimuon data together. Maximal radial flow is 
reached by the /?, due to its maximal coupling to pions, 
while all other hadrons follow some hierarchy in earlier 
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freeze-out. The T e ff values of the dimuon excess rise nearly 
linearly with mass up to the pole position of the p, but stay 
always well below the p line, completely consistent with 
the expectations for radial flow of an in-medium hadron- 
like source (here 7r + 7r~ — > p) decaying continuously into 
lepton pairs. 

This picture can still be refined. By modeling a p with 
the proper spectral shape and the mj spectrum as mea- 
sured, its contribution can be subtracted from the total 
measured distribution in the full M-m,T plane with the 
same side-window method as used for Fig. [9] below and 
described in [15]. The difference, essentially a continuum, 
is then refit, resulting in T c ff values for the "pure" in- 
medium (continuum) part in the 2w region. The modified 
plot is shown in the right panel of Fig. [5J The appearance 
is striking: all values are lower, but mostly so in the bin 
associated with the p pole (by ~20 MeV). This makes the 
shape even more sawtooth-like than before, and within 
errors the rise continues now up to about 1 GeV. 

Beyond the 2ir region, the T c g values of the excess 
dimuons show a sudden decline by about 50 MeV down 
to the IMR values. This decline is even more abrupt in 
the right panel of Fig. [8J than in the left and obviously 
connected to the in-medium emission itself, not to any pe- 
culiarities associated with the p peak. Extrapolating the 
lower-mass trend set by a hadron-like source to beyond 
1 GeV, such a fast transition is extremely hard to reconcile 
with emission sources which continue to be of dominantly 
hadronic origin in this region. A much more natural ex- 
planation would be a transition to a dominantly early, i.e. 
partonic emission source with processes like qq — > p + p~ , 
for which flow has not yet built up [L9ll23|. In this sense, 



the present analysis may well represent the first data- 
based evidence for thermal radiation of partonic origin 
in nuclear collisions, overcoming parton-hadron duality in 
the yield description on the basis of M-pT correlations as 
discussed in section 3. 

Theoretically, the extension of the unified LMR and 
IMR results over the complete M-px plane places severe 
constraints on the dynamical trajectories of the fireball 
evolution. Indeed all present scenarios [18, 19, 23, 26, 27] do 
not any longer rely on parton-hadron duality in the rates 
as in [21] . but explicitly differentiate between hadronic 
(mostly 47r) and partonic contributions in the IMR as al- 
ready discussed in connection with Fig. [3] The partonic 
fraction presently ranges from 0.65 for [26 (option EoS- 
B+ as used in Fig. [4} to "dominant" in [E ^ IST] . The 
exponential shape of the experimental tut spectra is re- 
produced by the models, consistent with the expectations 
for thermal radiation. However, due to remaining uncer- 
tainties in the equation-of-state, in the fireball evolution 
and in the role of hard processes [26], a quantitative de- 
scription of the much more sensitive mr-derivative T e ff 
vs. M in Fig. [5J is only slowly emerging. In particular, 
the more recent results from the authors of [19,23]l27]. 
while very encouraging in the description of the down- 
ward jump, are still preliminary and have not yet been 
formally published in their final form. A systematic com- 
parison of several model results to the data in Fig. [8J is 
therefore presently not possible. 
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Fig. 9. Excess yield ratios for peak, continuum and total vs. 
centrality for the mass window 0.2<M<1 GeV. Open charm 
is subtracted throughout. No acceptance correction applied. 
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tinuum, peak and total in the mass window 0.6<M<0.9 GeV 
(see also [4] right). Open charm is subtracted throughout. 



5 Centrality dependencies 

All data presented so far correspond either to (nearly) 
minimum bias, selecting dN^/ 'd?7>30, or to the semicen- 
tral window, selecting H0<dN c h/ dn<170; the results for 
the two conditions are very close. However, an enormous 
amount of information exists on the centrality dependence 
of practically every variable discussed in this paper. We 
select two topics of particular relevance here. 

The first one concerns the evolution of the shape of the 
excess mass spectra, following Figs. [2] and [3] In [15] . we 
have used both an rms analysis and a more sensitive side- 
window method to determine separately the yields of the 
peak and the underlying continuum. The centrality depen- 
dence from the latter is shown in Fig. [9] peak, continuum 
and total excess yield in the mass interval 0.2<M<1.0 
GeV, all normalized to the (fictitious) cocktail p with the 
assumption p/u=l (like in pp). The ui itself is directly 
measured, and its yield is found to be proportional to 
dN c h/dr]. The continuum and the total show a very strong 
increase, starting already in the peripheral region, while 
the peak slowly decreases from >1 to <1. Recalling that 
Fig. [9] is based on the excess mass spectra before accep- 
tance correction like in Fig. [3] roughly representing the 
full p spectral function, the excess/p ratio can directly be 
interpreted as the number of p generations created by for- 
mation and decay during the fireball evolution, including 
freeze-out: the "p clock", frequently discussed in the past. 
It reaches up to about 6 generations for central In-ln col- 
lisions; selecting low px this number doubles. 

The second topic concerns the centrality dependence 
of the slope parameter T e g for the excess data in the p-like 
window 0.6<Af<0.9 GeV, following Fig. H (right). Based 
as before on the side- window method [TSJ , the results are 
shown in Fig. [TO] separately for the p peak, the contin- 
uum and the total excess. The peak is seen to show a 
very strong rise, with hardly any saturation. However, the 
errors in the more central data become quite large, reflect- 
ing the continuously decreasing peak/total ratio as visible 



in Fig. [51 Conversely, continuum and total yield saturate 
much earlier. 

The large gap in T g between the peak and the con- 
tinuum, seen already in Fig. [7j has a much-debated in- 
teresting physics origin. The px spectrum of a thermal 
in- medium source is softer by a Lorentz factor M/E=l/~/ 
compared to that of a freely decaying freeze-out p [26 . On 
top, the in-medium values of T c g represent a temperature- 
flow average, while the freeze-out p receives the maximal 
flow. These two effects contribute about equally to the to- 
tal. The size of the gap ultimately reaches 70-100 MeV, 
but closes towards peripheral collisions. The w, with the 
same mass as the p, also shows a large gap to the p (com- 
pare Fig.rn}, which also closes finally to zero for the lowest 
pp-like window [IB"] . 



6 Evidence for cu in-medium modifications 

While most of the historical discussion on light-flavor vec- 
tor mesons in hot and dense matter has concentrated on 
the short-lived p (cr =1.3 fm), the longer-lived oj (23 fm) 
and <j) (46 fm) have received much less attention, since 
most of their dilepton decays occur after thermal freeze- 
out. Within the NA50 LMR analysis, the w and 4> have in- 
deed consistently been treated as "cocktail" particles and 
subtracted from the total (compare Fig. [I]). However, in- 
medium effects are expected for the (small) decay frac- 
tion inside the fireball, and these are actually contained 
in the Hees/Rapp scenario (see [26| and earlier references 
in there). 

NA60 has addressed the ui in a way directly coupled to 
the cocktail subtraction procedure. Due to the high mass 
resolution, the disappearance of the yield at low px out 
of the narrow to peak in the nominal pole position can 
sensitively be detected. The appearance of the yield else- 
where in the mass spectrum, originating from a mass shift, 
or broadening or both, is practically unmeasurable, due 
to the masking of the whole region by the much stronger 
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Fig. 11. Acceptance-corrected transverse mass spectra of the 
u> and the <f> for the semicentral bin. The T e g values correspond 
to different fit ranges for the two particles. A depletion of the uj 
yield at low mr relative to the fit line can clearly be recognized. 



7T7T — > p process, regenerating the p. Sensitive experiments 
on u> in-medium effects with clear clues as to their char- 
acteristics can therefore only be done in cold nuclear mat- 
ter experiments |29] . where p/uj = 1, but not in ultra- 
relativistic nuclear collisions. The evidence for the disap- 
pearance of cj's in the low tut region is shown in Fig. 1111 

As already mentioned in section 2, the u> and <j> are 
obtained as a byproduct of the cocktail subtraction proce- 
dure. The data are fit with the usual Tot exponential used 
before in Figs. |4] and El With respect to this reference line, 
there is hardly any anomaly visible for the </), but quite 
some loss for 1ow-pt oj's. The loss can be quantified with 
respect to the reference line, extrapolating down to zero. 
The fit parameters T c g for two different fit regions both for 
the uj and the (f> show the definition of the reference line 
to be quite uncritical. Forming the ratio data/reference 
line takes care of that part of radial flow which does not 
seriously affect the exponential slope. 

The results are shown in Fig. absolutely normal- 
ized as the full phase space ratio uj/N part . The effects of 
ui disappearance are quite striking: (i) a suppression of 
the relative yield below the reference line only occurs for 
Pt < 1 GeV; (ii) there is a very strong centrality depen- 
dence of the suppression, reaching down to < 0.5 of the 
reference line (the errors become huge for the central win- 
dow, because the uj can then hardly be recognized on top 
of the 7T7r processes at low my); (iii) the suppression effects 
are much larger than expected for the spectral distortions 
due to the blue shift from radial flow at low tot; a sim- 
ulation in the basis of the blast wave parameters from 
[2] shows at most 10% effects for central collisions. The- 
oretical simulations addressing these results are not yet 
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Fig. 12. pt dependence of the u yield with respect to the fit 
line in Fig. 1111 absolutely normalized for the full phase space 
yield, for different centralities. The solid line for pr < 1 GeV 
shows the result from a blast-wave fit to the uj for central col- 
lisions. The dotted lines are only meant to guide the eye. The 
errors are purely statistical. The systematic errors are negligi- 
bly small compared to the statistical ones. 



available. It should be added that the same procedure ap- 
plied to the 4> does describe the data solely on the basis 
of radial flow. No effect beyond that can be recognized, 
within errors. 



7 Conclusions 

This paper, supplementing [14] . contains the most com- 
prehensive data set on excess dileptons above the known 
sources which has so far become available through NA60. 
We have concentrated here more than before on interpre- 
tational aspects, in particular on the way, " parton-hadron 
duality" in the yields can be overcome by a careful study 
of M~pt correlations. The data mediate a clear conclusion 
on the dominance of partonic processes for M>1 GeV. 
A systematic comparison with theoretical models reveals 
remaining ambiguities in the modeling, but the overall 
agreement with the data tends by now to support the 
same conclusion. 
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